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The hypothesis of a discrete fabric of the universe-the “Planck scale”-is always on stage, since it 
solves mathematical and conceptual problems in the infinitely small. However, it clashes with special 
relativity, which is designed for the continuum. Here we show how the clash can be overcome within a 
discrete quantum theory where the evolution of fields is described by a quantum cellular automaton. 
The reconciliation is achieved by defining the change of observer as a change of representation of the 
dynamics, without any reference to space-time. We use the relativity principle, i.e. the invariance of 
dynamics under change of inertial observer, to identify a change of inertial frame with a symmetry of 
the dynamics. We consider the full group of such symmetries, and recover the usual Lorentz group 
in the relativistic regime of low energies, while at the Planck scale the covariance is nonlinearly 
distorted. 


I. INTRODUCTION 

Is the world continuous or discrete? Richard Feynman 
m® motivated a discrete universe as the only way it can 
be simulated by its own constituents, which means by a 
quantum computer. Einstein himself considered a dis¬ 
crete space-time as a possibility, however, he complained 
about the lack of an appropriate mathematical frame¬ 
work 0. Usually we dismiss discreteness on the basis of 
a mathematical convenience of continuos theories. But 
the continuum leads to still unsolved mathematical prob¬ 
lems in the infinitely small, problems that do not arise in 
the discrete. The discrete, on the other hand, seems to 
raise a problem: the disagreement with Einstein’s special 
relativity. 

The debate about the clash between a discrete space 
and Lorentz symmetry has been recently renewed be¬ 
cause in some approaches to quantum gravity (such as 
Regge Calculus [4], spin-foam [5], causal sets, 0 ) the 
fundamental description of space-time is a discrete struc¬ 
ture, to which the continuum is only an approximation. 
The scale of this discreteness is the Planck length, which 
is amazingly small-the Planck length compared to a me¬ 
ter is like the electron radius compared to the size of our 
galaxy. 

Why the Lorentz transformations would not work with 
a discrete space-time? The objection is that a discrete 
space-time would not be invariant under the Lorentz 
group, even if we take it as discrete. Such a point, raised 
more than sixty years ago, was disproved for d = 3 space 
dimensions in Ref. [7 : . However, it was shown that the 
minimum admissible boost would be huge: 0.866 times 
the speed of light! It seems therefore that there is no 
way for the reconciliation of Einstein’s special relativ¬ 
ity with a discrete fabric of space-time. The issue is, 
however, a false problem, originating from the unneces- 
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sary requirement of enforcing a covariance designed for 
the continuum. The right point of view is to take the 
Lorentz covariance only as an approximate symmetry, 
and recovering it in the regime where the discreteness 
looks continuous. This is similar to what happens for a 
crystalline medium, which looks isotropic at large scales, 
whereas instead it is highly anisotropic at microscopic 
scales. The smaller is the crystal structure, the more ac¬ 
curate is the continuum symmetry: think that the Planck 
length is 10 -25 A! 

The right point of view is thus to consider the contin¬ 
uum as an approximation of the discrete when observed 
at very large scales. It is thus conceptually legitimate 
that the Lorentz transformations are actually distorted 
at the tiny Planck scale. An example of such a distorted 
Lorentz symmetry is that of doubly special relativity it 
uni, where the distorted Lorentz transformations, in ad¬ 
dition to the speed of light, preserve also an energy scale. 

Here we show that if we take only the very essence of 
the relativity principle-the invariance of the physical law 
under change of inertial representation-we get nonlin¬ 
ear Lorentz transformations, which happen to be of the 
same kind as those of doubly special relativity. In the 
continuum description the reference frame is a Cartesian 
coordinate system-what is called “position representa¬ 
tion” in quantum theory. Other representations of the 
dynamics are given in terms of constants of motion-such 
as the momentum or energy represent at ions-and these 
provide a viable notion of inertial frame in a world made 
of countably many quantum systems. 

II. RESULTS 

We take the pragmatic point of view that quantum 
theory is more efficient a description of our world than 
classical theory, pretty much like Kepler laws for the 
planet orbits are more efficient than Ptolemaic epicycles. 
We therefore consider the most general quantum discrete 
theory, which is a quantum cellular automaton mmm- 
This consists of infinitely many quantum systems (qubits, 
Fermionic or Bosonic fields), whose evolution occurs in 
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FIG. 1: The green surface represents the dispersion relation 
in Eq. ^ where we fixed k z — 0 (c^k = arccos (k) have the 
same plot). The red surface is the usual relativistic dispersion 
relation = k x + ky. Notice that the two surfaces get closer 
approaching the origin for k —> 0. 

discrete steps, and which interact locally , namely every 
system interacts with a bounded number of other sys¬ 
tems. A consequence of locality is that signals propagate 
at finite speed over the interacting network. For our pur¬ 
poses it is sufficient to consider a single particle, and 
technically this simplifies the automaton to be linear in 
the quantum field. Moreover we require that the dynam¬ 
ics is reversible, hence it is described by a unitary matrix: 
this is what is called Quantum Walk [13, 14]. 

In the discrete context universality of the physical law 
corresponds to the homogeneity of the evolution. This 
implies that the network of interactions represents a 
group G of “translations”, moving from a system to an¬ 
other interacting with it: this is a so-called Cayley graph 
of G (see Methods). Moreover, if we want to reproduce 
the physical situation of flat space where parallel trans¬ 
ports along two different paths end up in the same result, 
then the “translation” group must be Abelian. This al¬ 
lows us to make use the Fourier transform. Upon denot¬ 
ing by |k) the eigenvector of the translations correspond¬ 
ing to wave-vector k, we write the quantum walk unitary 
operator A as follows 

A = j dk^ k <g> |k)(k|, (1) 

where B denotes the Brillouin zone of the interaction 
lattice, and Ak is a unitary 5 dimensional matrix, with s 
finite. As proved in Ref. m, the easiest nontrivial quan¬ 
tum walk has s = 2. The eigenvalues of A are obtained 
by solving the eigenvalue equation 

A k tp(u, k) = e lu ip(w, k) (2) 

which can be rewritten in relativistic notation as follows 

n fl (k)a fJ 'i/;(k) = 0, (3) 

where we introduced the four-vectors k = (o;,k), n(k) = 
(sincj, n(k)), and cr = (/, cr), with cr = (<r x , cr y , a z ) de¬ 
noting the Pauli matrices, and the vector n(k) are defined 


by 

n (k) • cr ■= - -4c)- ( 4 ) 

The eigenvalues can be collected into two functions 
cj ± (k) called dispersion relations. In this scenario the 
constants of motions are k and hence a change of 
representation corresponds to a map k i-T fc'(fc). Now 
the principle of relativity corresponds to the requirement 
that the eigenvalue equation § is preserved under a 
change of representation as follows 

n M (fcK = f^n M (fcVr fc , (5) 

where T^, T& are invertible matrices. 

Eq. (J5| translates the relativity principle for the QW 
evolution: the dynamics is left invariant by a change of 
observer. 

The simplest example of change of observer is the one 
given by the trivial relabeling k' = k and by the matri¬ 
ces T/c = Tk = e zA ( k ), where A(k) is an arbitrary real 
function of k. When A(k) is a linear function we recover 
the usual group of translations. The set of changes of 
representation k \-> k'{k) for which Eq. © holds are a 
group, which is the largest group of symmetries of the 
dynamics. 

If to the general assumptions defining the quantum 
walk we just add that of isotropy, it turns out that there 
are only two admissible quantum walks HE which in 
the small wave-vector regime give exactly the two Weyl 
equations for the left and right massless Fermion. Indeed, 
with the above assumptions the only possible lattice is 
the body centered cubic one, and modulo local unitary 
equivalence the two admissible quantum walks are 

A ^ := A ± (k )I — m ± (k) • cr ± , (6) 

where 

( S X CyC Z A C X SyS Z \ 

c x s v c z =F S x c y s z , 

CxCy^z A S x SyC z J 

X (k) := ( c x CyC z T s x SyS z ), (7) 

c a := cos (k a /y/3), s a := sin(/c a /v / 3), a = x,y,z, 

where <r + = cr and a = a r , with T denoting the trans¬ 
posed matrix. The dispersion relations are given by 

n^(fc)n /i± (fc) = 0, (8) 

and are plotted in Fig. [1] 

In the small wave-vector regime k ~ ko = (0, 0, 0) one 
has n(k) ~ k , recovering the usual relativistic disper¬ 
sion relation. The Weyl equations can be also recovered 
in the neighborhood of the wavevectors ki = |(1,1,1), 
k 2 = —1(1,1,1), k 3 = —1(1,0, 0). The mapping be¬ 
tween the vectors k^ exchange chirality of the particle 
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FIG. 2: The distortion effects of the Lorentz group for the discrete Planck-scale theory represented by the quantum walk in 
Eq. ©• Left figure: the orbit of the wavevectors k = (k x , 0, 0), with k x G {.05, .2, .5,1,1.7} under the rotation around the z 
axis. Right figure: the orbit of wavevectors with |k| = 0.01 for various directions in the (k x ,k y ) plane under the boosts with (3 
parallel to k and \/3\ G [0,tanh4]. 


and double the particles to four species in total. There¬ 
fore we have four different particles-two left-handed and 
two right-handed-namely the discreteness also doubles 
the particles, which is the well known phenomenon of 
Fermion doubling H5|. In the following the term “small 
wavevector” will denote the neighborhoods of the vectors 
k; i = 0,... 3. 

We now show that the group of symmetries of the dy¬ 
namics of the quantum walks in Eq. © contains a non¬ 
linear representation of the Poincare group, which ex¬ 
actly recovers the usual linear one in the small wave- 
vector regime. For any arbitrary non vanishing function 
f(k) we can introduce the four-vector 

p(f) _ x>(f\k) := /(fc)n(fc) (9) 

and rewrite the eigenvalue equation © as follows 

p(f } <j fJ '!p(k) = 0. (10) 

Upon denoting the usual Lorentz transformation by Lp 
for a suitable / (an example is provided in the supplemen¬ 
tal material) the Brillouin zone splits into four regions 
i = 1,..., 4 centered around k^ i = 0,... 3, such that the 
composition 

C { p := ( 11 ) 

is well defined on each region separately (see Methods). 
The four invariant regions corresponding to the four dif¬ 
ferent massless Fermionic particles show that the Wigner 
notion of ” particle” as invariant of the Poincare group 
survives in a discrete world, consistent with a physical in¬ 
terpretation of the Fermion-doubled particles. For fixed 
function / the maps provide a non-linear represen¬ 
tation of the Lorentz group mmm- In Figs. [2] and 
Owe show the numerical evaluation of some wavevector 


FIG. 3: The green surface represents the orbit of the wavevec¬ 
tor k = (0.3, 0,0) under the full rotation group 50(3). 


orbits under subgroups of the nonlinear Lorentz. The 
distortion effects due to underlying discreteness are ev¬ 
ident at large wavevectors and boosts. The relabeling 
k fc'(fc) = £p\k) satisfies ^ with T^ = Ap and 
Tk = Ap for the right-handed particles, and T/- = Ap and 
Tk = Ap for the left-handed particles, with Ap and Ap 
being the (0, |) and (^,0) representation of the Lorentz 
group, independently on k in each pertaining region. 

For varying /, we obtain a much larger group, includ¬ 
ing infinitely many copies of the nonlinear Lorentz one. 
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In the small wave-vector regime the whole group collapses 
to the usual linear Lorentz group for each particle. 

Up to now we have analyzed what happens with mass¬ 
less particles. A simple way to obtain the Dirac equation 
is to pair an automaton in Eq. ([6| with its adjoint into a 
direct sum, as in Ref. HE thus leading to a new automa¬ 
ton giving the Dirac equation in the small wave-vector 
regime. A relevant feature of the discreteness is that be¬ 
cause of unitarity the mass parameter is upper bounded 
[XH1 • Now if one derives the full symmetry group of the 
dynamics as we have done for the two automata in Eq. 
© one discovers that the group is a nonlinear representa¬ 
tion of the de Sitter group 50(1,4) with infinite cosmo¬ 
logical constant, with the rest mass of the particle playing 
the role of the additional coordinate (see Methods). It is 
noticeable that even for pure boosts the rest-mass is in¬ 
volved in the transformation. For rest-mass much smaller 
than the upper bound and for pure boosts one recovers 
the previous nonlinear Lorentz group for zero-mass. 

III. CONCLUSION 

We have seen what happens of the Lorentz group in 
a quantum world that is discrete. The main point is to 
abandon the idea of enforcing the exact Lorentz symme¬ 
try on the discrete, but instead to consider the symmetry 
as an approximate one that holds only in the small wave- 
vector and small mass regime. But the natural question 
is now: how small? According to the common opinion the 
scale of discreteness a is identified with the Planck scale. 
In terms of the maximum wavevector kM in the Bril- 
louin zone, one has In the small wavevector 

regime we recover the simple relations [181 c = and 
h = /iac, with c, /i, and r denoting the speed of light, the 
Planck constant, and the time-step, respectively. Then 
the maximum mass /i of the quantum walk is the Planck- 
mass. A way of deriving /a and a heuristically is to keep 
literally the argument of taking the mass of the parti¬ 
cle bounded in order to keep the Compton wave-length 
A c larger than the Schwartschild radius. Noticeably for 
m = /I the dispersion relation is constant, namely with 
no propagation of information, a situation reminiscent of 
a micro black hole m ■ Remarkably general relativity en¬ 


ters the present quantum digital framework also through 
the unforeseeable appearance of the De Sitter symmetry 
group, which connects different Dirac particle mass val¬ 
ues. Are these only coincidences? The dream is that it 
is a new route to quantum gravity. 

The crucial question is now what can be actually seen 
experimentally. Recently, experimental tests of Planck- 
scale phenomenology have been proposed [201422] . In par¬ 
ticular, the modification to the usual dispersion relations 
can in principle be detected in observation of gamma-ray 
bursts from deep-space events, where billions of light-year 
of distance can sufficiently amplify the weak vacuum dis¬ 
persive behavior due to discreteness [23]. In our context 
this can be proved with the free electromagnetic field de¬ 
rived as the two-particle sector of the quantum walk in 
Eq. ([6| [24]. This possibility reconnects with the recent 
analysis of data [25] from Fermi-LAT concluding that the 
observations set an upper-bound to the scale of discrete¬ 
ness which is smaller than the Planck scale a by a factor 
2.8. The analysis of Ref. [25] can be refined with a com¬ 
plete theoretical derivation based on Ref. [24] and on the 
results presented in this letter. This would also take into 
account the possibility of a compensating effect due to 
the phenomenon of relative locality [26]. In short rela¬ 
tive locality is the phenomenon due to the nonlinearity 
of the Lorentz transformations, which generalizes the rel¬ 
ativity of simultaneity to relativity of the full space-time 
coincidence of events. The separation of events under 
boost is amplified by the difference of their frequency 
domain. Indeed the Fermi-LAT observation is based on 
a predicted time-delay between two events with a huge 
difference in frequency, which could then be compensated 
by the relative-locality effect. The fully fledged discrete 
theory given here, derived from very general principles, 
allows for a thorough quantitative evaluation that takes 
into account both the dispersive vacuum and nonlinear 
Lorentz transformations. 
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SUPPLEMENTAL INFORMATION 


IV. QUANTUM WALK FROM PRINCIPLES 

Quantum Walks (QWs) describe the evolution of a 
quantum particle over a lattice. The dynamics is as¬ 
sumed to be reversible , hence the QW will be represented 
by a unitary transformation. By denoting with G the set 
of the lattice’s points we can convenienty introduce the 
Hilbert space I 2 {G) and the orthonormal basis |x) which 
corresponds to the position of the particle. If we associate 
the Hilbert space C s to the internal degrees of freedom 
of the particle, the QW is then a unitary operator on 
£ 2 (G) ® C s . 

a. Locality The QW evolution is assumed to be lo¬ 
cal i.e. information propagates through the lattice at a 
bounded speed. Given a lattice, let 7V X be the set of 
nearest neighbors of the site x. If the particle is localized 
at site x, then, after one step of the QW, it must be lo¬ 
calized within a finite set 7V X . Such a locality condition 
introduces a notion of causal cone in the lattice. 

b. Homogeneity The evolution is assumed to be ho¬ 
mogeneous , i.e. requirement that all the sites are equiv¬ 
alent. The QW evolution should not not allow one to 
discriminate between two points x and x'. This require¬ 
ment (see Ref. 1 for a full derivation and Appendix |A| 
for a shorter review) implies that the set og points G is a 
group and the lattice is a possible Cayley Graph of this 
group. 

c. Isotropy The assumption of isotropy translate 
the requirement that there is no priviledged direction on 
the lattice. The mathematical traslation of this require¬ 
ment requires the existence of a group of permutation 
that act on the generators of the group G that can be 
faithfully represented on the internal degrees of freedom 
(see Ref. |T] for a full derivation and the Appendix |A| for 
a shorter review) 

d. Flat and curved space The above sketched frame¬ 
work encompass a broad variety of dynamics. In partic¬ 
ular, depending on the properties of the group G, we can 
have a quantum dynamics on a generally curved space. 
If we are interested to make contact with special relativ¬ 
ity, it is natural to restrict the scenario to QWs corre¬ 
sponding to dynamics on the flat Minkovski spacetime. 
This requirement correspond to assume the group G to 
be virtually abelian , i.e. G has an abelian subgroup of 
finite order. We can further restrict ourselves to the case 
in which G is abelian without any loss of generality. The 
price to pay for this restriction is to add additional inter¬ 
nal degrees of freedom (see Ref. [2J and Appendix |A| for 
a more complete discussion). 

e. Fourier analysis If the group G is abelian it is 
convenient to study the dynamics in the Fourier trans¬ 
form basis |k) := (2tt) 2 e zk ' x |x). Since homogene¬ 
ity condition implies that the QW commutes with the 
translations on the lattice, in the Fourier basis the QW 


operator can be written as 

A = f d 3 k |k)(k| G A k ( 12 ) 

Jb 

where B denotes the first Brillouin zone of the underlying 
lattice. The unitary constraint implies that Ak is unitary 
for every k G B and the locality assumption implies that 
Ak is a matrix polynomial in e zh k . Notice that due to the 
discreteness of the lattice the QW is band-limited in k. 
The spectrum {e tuJ ^ } of the operator Ak, and expecially 
its dispersion relation that is the expression of the phases 
as functions of k, plays a crucial role in the analysis 
of the QW dynamics. 

/. Weyl Quantum Walk If the dimension of the 
Hilbert space of the internal degrees of freedom is s = 2 
and the group G is Z 3 the requirements of locality , ho¬ 
mogeneity and isotropy then [I] the QWs can only be 
defined over the body-centered-cubic lattice and they are 
equivalent (up to a local change of basis) to the following 
two QWs: 

A^ := A ± (k)/ — m ± (k) • cr ± , (13) 

( Sx^yCz d: C x SyS z \ 

Cx^yCz “F SxCy^z I 5 

CxCy^z d: S x SyC z J 

A (k) := ( c x CyC z G- s x SyS z ), (Id) 

c a := cos(& a /\/3), s a := sin(fc a /\/3), a = x,y,z. 

The Pauli matrices cr + = a are the usual ones, while the 
cr~ = cr T are their transposed ones, and k G B where B 
denotes the Brillouin zone of the BCC lattice. 


V. NONLINEAR LORENTZ 
TRANSFORMATIONS FOR THE WEYL 
QUANTUM WALK 

In this section we give the explicit construction of the 
non-linear representation of the Lorentz group on the set 
of the solutions of the Weyl QW dynamics. The set of 
solution will be split into four regions, each one of them 
carrying a non-linear deformation of the Lorentz group. 
Let us consider the splitting of the Brillouin zone B 

Bq := {k G B|A(k) > 0, cos (2 k y /Vs) > 0}, 

B[ := {k G B|A(k) < 0, cos (2 k y /V?>) > 0}, 

/— (15) 

B' 2 := {k G B|A(k) > 0,cos (2k y /Vs) < 0}, 

B 3 := {k G B|A(k) < 0, cos (2 k y /Vs) < 0}, 

with B = U 3 = 0 B' up to a nullmeasure set, and let us 
denote with nW(k) the restriction of n(k) to B'. Notice 
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that we dropped the label gfc denoting the two different 
Weyl walk, since the results holds the same in the two 
cases. We now denote with U the unit open ball in M 3 (U 
denotes its closure) and with H the subset of U defined 
as follows 

H := {m G U s.t.m x = ±m z , 2 m 2 + 2ra 2 > 1}. (16) 

We then consider the regions 

B< -nW-^UXH) (17) 

and the function /(u;, k) defined as follows: 

/(cj,k) = 5 f (n(k)) = 

=g(r,0,(j)) := 1 + r[ds (+ 1 ) 

a(r) := 1 — r 2 , 

b(r,0,</>) := (cos 2 ( 0 ) -sin 2 ( 0 )) 2 + 

+ (| — r 2 (l - cos 2 ($) sin 2 ( 0))) 2 

(18) 

where we used the spherical coordinates n x ( k) = 
r cos# cos 0, n y ( k) = rsin#, n z ( k) = rcos#sin(/>. Fi¬ 
nally, we define the maps as 

: -> r ». i,<<>: (k) ~ /("• k > („ s (<" ( k)) ■ 

Z* := {(a;, k) s.t. k G Eh, sin 2 c o — |k| 2 = 0}, 
r 0 := {p G M 4 s.t. p^ = 0}. 

(19) 

One can prove (see Appendix ) that the maps define 
an analytic diffeomorphism between Z^ and To. As a 
consequence, the composition 

cf : Z; Z; cf := V- 1 oLpoV (20) 

is a well defined nonlinear representation of the Lorentz 
group on the set Z im Since the union of the sets coin¬ 
cides with the whole (up to a null measure set) Brillouin 
zone, we have that the collection of the maps provide 
a notion of Lorentz transformation for any (up to a null 
measure set) solution of the Weyl QW dynamics. 

VI. DYNAMICAL SYMMETRIES OF DIRAC 
QUANTUM WALK 

The Dirac quantum walk is obtained by performing 
the direct-sum of two Weyl walks and introducing off- 
diagonal blocks in such a way that the obtained matrix 
is a well defined quantum walk. As proved in Ref. PQ, 
there are only two admissible walk 

D± _ (nA ± iml \ 

_ yiml nA ± J (21) 

0 < n, m < 1 n 2 + m 2 = 1 


The eigenvalue equation of the Dirac QW can be written 
as 

\pi f) (w, k, m)Y - ml]ip(u, k, m) = 0 , ( 22 ) 

where are the Dirac 7 matrices in the Weyl representa¬ 
tion, and m is then interpreted as the particle mass. Due 
to the explicit dependence of from m the covariance 
under change of reference cannot leave the value of m 
invariant. In such case the dispersion relation resorts to 
the conservation of the de Sitter norm 

sin 2 uj — (1 — m 2 ) |n(k) | 2 — m 2 = 0. (23) 

The group leaving Eq. ( [23] ) invariant is the De Sitter 
group 50(1,4). In the limit of m <C 1, the usual Lorentz 
symmetry is recovered. The analysis of De Sitter covari¬ 
ance of Eq. ( [22] ) will be given in a forthcoming publica¬ 
tion. 


Appendix A: 3-dimensional Quantum Walks and the 
Weyl Quantum Walk 

We consider Quantum Walks (QWs) that describe the 
evolution of a quantum particle over a lattice of M 3 . The 
dynamics is assumed to be reversible, hence the QW will 
be represented by a unitary transformation. We consider 
Bravais lattices , i.e. lattices that are generated by linear 
independent discrete translations. By denoting with G 
the set of the lattice’s points we can convenienty intro¬ 
duce the Hilbert space £ 2 (G) and the orthonormal basis 
|x) which corresponds to the position of the particle. The 
internal degrees of freedom of the particle are described 
by the Hilbert space C s . 

The QW evolution is assumed to be local i.e. informa¬ 
tion propagates through the lattice at a bounded speed. 
Given a Bravais lattice, let N x be the set of nearest neigh¬ 
bors of the site x. If the particle is localized at site 
x, then, after one step of the QW, it must be localized 
within the set 7V X . The set 7V X corresponds to set of vec¬ 
tors 5 := {h} such that 7V X = {y s.t. y = x + h, h G 5} 
Such a locality condition introduces a notion of causal 
cone in the lattice. 

The QW evolution is also assumed to be homogeneous , 
i.e. it must commute with the discrete traslations. Let 
Xh be the translation operators on ^ 2 (G), such that i.e. 
Th|x) = |x + h) These assumption implies that the QW 
evolution can be represented by a unitary operator on 
£ 2 (G) 0 C s as follows: 

A:=Y,T h ®A h (Al) 

h es 

where A h are s x s matrices and Th is the translations 
that connect the site x with its nearest neighbor x + h. 

The QW dynamics is conveniently study in the Fourier 
transform basis |k) := ( 2tt )a X] x e2k x | x )’ which the 



QW operator A becomes 

A = f d 3 k |k)(k| <S> A^, A k := ^ e~ tkh A h , 

hes 

(A 2 ) 

where B denotes the first Brillouin zone of the underly¬ 
ing lattice and it is defined by the following linear con¬ 
straints: 

B := P|{k G R 3 | - 7r|hj| 2 < |k • h*| < TrlyJ 2 } 

i (A3) 

h i ' hj = Sij , h j G S. 

The unitary constraint implies that A^ is unitary for ev¬ 
ery k G B and the locality assumption implies that A^ 
is a matrix polynomial in e zhk . Notice that due to the 
discreteness of the lattice the QW is band-limited in k. 
The spectrum {e 2U T } of the operator Ak, and expecially 
its dispersion relation that is the expression of the phases 
as functions of k, plays a crucial role in the analysis 
of the QW dynamics. 

It was shown in Ref. jli that, when the dimension of the 
Hilbert space of the internal degrees of freedom is s = 2 , 
the additional assumption of isotropy A of the evolution 
implies that the QWs can only be defined over the body- 
centered-cubic lattice and they are equivalent (up to a 
local change of basis) to the following two QWs: 

:= A ± (k)J — in ± (k) • a ± , (A4) 

where we define 

( SxCyCz i CxSySz\ 

C x SyC z T SxCy^z I i 
^x^ySz A: S x SyC z J 

A ± (k) := (c x c y c z =F s x s y s z ), (A5) 

c a := cos(fc Q /v / 3), s a := sin(fc a /v / 3), a = x,y,z. 

The Pauli matrices cr+ = a are the usual ones, while the 
<j~ = <j t are their transposed ones, and k G B where B 
denotes the Brillouin zone of the BCC lattice (see Fig. [4]). 
In the limit |k| <C 1 we have 

n± ( k ) ~ A i ~ ex p[ _i ^ ■ cr± ]’ ( A6 ) 

corresponding to the evolution of the Weyl’s equation 
with the rescaling k. Since the A + and A~ re¬ 

produce the dynamics of the Weyl equation in the limit 
|k| < 1 , we refer to them as Weyl QW. 

For the sake of simplicity, in the following we will con¬ 
sider only one Weyl QW, i.e. we define A^ := Aj, the 
choice of Pauli matrices is the usual one cr := cr + , and 
we similarly drop all the db superscripts. This choice is 
completely painless since all the results of this paper can 
be easily adapted to the choice Ak = A^. In order to 
simplify the notation we also adopt the convention 

*■ ^ % (A7) 



FIG. 4: Top : Body centered cubic lattice. Each pair of near¬ 
est neighbor is connected by a blue link. Each point of the 
lattice has 8 nearest neighbors, e.g. the red points are the 
nearest neighbors of the yellow one. Bottom : Brillouin zone 
B of the BCC lattice. The zone is a rhombic dodecahedron in 
which the opposite faces are identified. 


in order to get rid of the annoying ^ factor. With this 
choice we have that in the |k| <C 1 limit, n(k) ~ k. 


Appendix B: Deformed relativity and Weyl 
quantum walk 

In this section we prove that the construction in Eq. 
is a well defined deformed Lorentz symmetry for each set 
YLi Let us define 

V = M o V, V:(u>, k) i-> (w, n(k)) (Bl) 

where we also assumed /(cj, k) = g(u), n(k)). 

We now need to study separately the properties of the 
two maps V and J\f. 
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1. Study of the map n(k) 


In this subsection we study the anal ytic al properties of 
the map V, which, according to Eq. (Bl) resorts to the 
map n. The analysis will proceed through the determina¬ 
tion of the largest subdomains {B-} of invertibility of the 


map n. We first prove that on the closure of each do¬ 
main the map is surjective on the closed unit sphere U. 
Then, we determine the geometry of the ranges n(B'), 
showing that they are homotopic to a solid genus-two 
torus. 

Let us denote by B the Brillouin zone of the center cu¬ 
bic lattice. B, upon a proper identification of its bound¬ 
ary points (see Fig. [4| is a compact 3-dimensional man¬ 
ifold. The Jacobian J n (k) of the map n(k) is given by 


J n (k) := det[< 9 * 7 ij(k)] = cos( 2 A: ?/ )A(k), (B 2 ) 

and it vanishes on the set 


F = G U X, 

X := {k E B| cos(2 k y ) = 0}, (B3) 

G := {k E B|A(k) = 0}. 

Since VA(k) fiz 0 for all k such that A(k) = 0, the implicit 
function theorem guarantees that G is a well defined 2 - 
dimensional surface. In the following we will denote by 
{B'} (* ranging in some set) the disjoint connected sub¬ 
sets of B \ F, thus 

B \ F = |J B', B' n B' = 0 for i ± j. (B4) 

i 

For each i the set B' is open and we denote as B' its 
closure and as dB[ its boundary. 

Now let us denote with U CM 3 the closed unit-radius 
ball, and with S 2 the sphere of radius 1 in M 3 . Moreover, 
let us define the parametric curves 

e±(t) : = -^=(sin(t), cos(t), ± sin(t)) T (B5) 
v 2 


For each i, let nW(k) denote the restriction of the map 
n(k) to the set B-. Then nW(k) defines a diffeomorphism 
between B^ and its image nW(B-) and we have 

n(°)(B')=n( 2 )(B') = Q 0 

n«(B' 1 )=n( 3 )(B') = Q 6 . 

The proof of this result is rather involved and can be 
found in Appendix |C| The B' regions are plotted in Fig. 
[5] The most important consequence of this result is that, 
for each i, the set n(B-) (see Fig. [ 6 | i) coincides with U 
except a null measure set and ii) it is homeomorphic to 
a genus two torus. 


2. Study of the map AT 

Since for all i the region n(B') has a nontrivial topol¬ 
ogy, the set {(cj,m) s.t. \uj\ < |,m E n(B'), sin 2 uj — 
|m | 2 = 0} cannot be diffeomorphic to any Lorentz- 
invariant region of M 4 . A possible way to change the 
topology of n(B') is to exclude the set H C U (as it is 
shown in Fig. pi) of vectors m satisfying the following 
inequalities 


' rn x = ±m z , 

< 2 m 2 + m 2 < 1 , (B9) 

„ 2m 2 + 2m 2 > 1. 

Then, the set U \ H is topologically trivial and we have 
U\H C n(B-) for all i. Let us now consider the function 
J\f : (u;,m) (p 0 ,p) = ^(cj, m)(sino;, m) restricted to 

the set 

N := {(cj, m) s.t m E U \ H, |cj| < sin 2 c o — |m | 2 = 0}. 

(B10) 

As shown in Appendix [D| it is possible to define the func¬ 
tion g(uj, m) is such that Af defines a diffeomorphism be¬ 
tween N and the null mass-shell 


and the sets 


Qa 

Qb 

Ti 

T 2 


U\(e + (T 1 )Ue_(T 2 )) 
U \ ( e +(T 2 ) U e_(T 1 )) 

/_ 7T 7T \ 

V 2 ’ 2 / 

(-7T,-f)U(f,7r]. 


(B6) 


Given all the definition introduced in this section, we 
have then the following result: 


r 0 := {p e M 4 , s.t. p'-'p/j = 0}. (Bll) 

and that its Jacobian matrix at tho origin is 0, i.e 

JaKO) = I. (B12) 

Finally, for each i we denote by B^ the counter-image 
of the set U \ H under the map nW and by the 
composition 


Lemma 1 There are four different B- regions, deter¬ 
mined by the following conditions 

Bq := {k E B|A(k) > 0, cos(2 k y ) > 0}, 

B i := {k G B|A(k) < 0, cos(2 k y ) > 0}, 

Bf, := {k G B|A(k) > 0, cos(2 k y ) < 0}, 1 ’ 

Bg := {k G B|A(k) < 0,cos(2fc y ) < 0}. 


V {i) : Zi -p r 0 V {i) := Af o V (i) (B13) 

P<‘> : I, -, N : (") « („ ( ,r (k) ) 

A/”: N —> To V :^) -> J/(w,m) ^ S1 ^, W ) 


Ti := {(cj, k) s.t. k E B^, sin 2 uj — |k | 2 = 0}. 
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FIG. 5: (Colors online). Top left figure: surfaces A(k) = 0 in Eq. ( |A5| ) (yellow) and cos(2 k y ) = 0 (red planes) inside the 
Brillouin zone (transparent). Top middle figure: B 0 zone (red X-shaped). Top right figure: B 0 (red) and Bi (blue). Bottom 
left to right: Bi, B2, B3. Bottom right: region Bi represented in a properly translated Brillouin zone. 

In this paper the Lorentz transformations are those that leave the dispersion relations of the Weyl QW invariant, and act on 
the Weyl spinor independently of the wave-vector. In such way they are nonlinear in (u, k) and linear over the Weyl spinor. 
Therefore the Lorentz group acts as a group of diffeomorphisms over the Brillouin zone B. The four domains Bi C B are Lorentz 
invariant (up to a null-measure set, see Fig. [6j . More precisely a point (cj, k) with k 6 Bi and sin 2 u — |n(k)| 2 = 0 is mapped 
to a point (o/jk') with sin 2 uj' — |n(k')| 2 = 0 and k x £ B*. Moreover, the map n maps each B* into the same set (up to null 
measure set: see Fig. [6j . Since the kinematics of a wavevector k depends only on the vector n(k), we can conclude that the 
Bi regions are kinematically equivalent and they can be interpreted as four different massless Weyl Fermions. Because of the 
identification of the boundary points in the Brillouin zone, all the Bi regions have the same X-shape as Bo- This is evident in 
the bottom right figure, in which we see that the region Bi (in red), when represented in a properly translated Brillouin zone 
(in blue), has the same X-shape as the region Bo. Considering the identification of the boundary points of the Brillouin zone 
in Fig. R] one realizes that the opposite arms of the X are glued together, resulting in a solid double-torus (genus-two). This 
result is rigorously proved in the text where we show that the B* regions are diffeomorphic to a solid ball pierced by two arches 
of ellipses (Fig. [6|. 


For each i, the map [ s an analitic diffeomorphism 
between the region X* and the Lorentz invariant set l” 0 
which satisfies the condition Jx> i ( 0 ) = I. Then the com¬ 
position 

4 ° : -»• Ei Cf := 2T 1 oLpoV (B14) 

is a well defined nonlinear representation of the Lorentz 
group on the set X^. Since the union of the B^ sets coin¬ 
cides with the whole (up to a null measure set) Brillouin 
zone, we have that the collection of the maps provide 
a notion of Lorentz transformation for any (up to a null 
measure set) solution of the Weyl QW dynamics. 


Appendix C: Proof of Lemma [l] 

In this section will give the proofs of the results con¬ 
tained in Lemma [l] Since the proof is quite involved, we 
split it into several pieces. Let us begin by defining the 
sets 

Q' := U \ R R := S 2 U E + U E_ . (Cl) 

Obviously Q' is open and connected, with Q' = U and 
dQ' = R. We now prove some useful properties of the 
map n. 
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FIG. 6: (Colors online) Left figure: region Q a . Right figure: H zone in red inside the unit ball. In the left figure, the tubes 
around the arches e+(T)i and e_(T)2 emphasize the piercing of the ball by the one-dimensional holes along the elliptic arches 
e+(T)i and e_(T)2. The region Q a is clearly homeomorphic to a solid toms of genus two. Because of this non-trivial topological 
feature the set {(ca,m) s.t. \uj\ < f ,m G n(B*), sin 2 u — |m| 2 = 0} cannot be diffeomorphic to any Lorentz-invariant region of 
M 4 . However it is possible to remove from the region Q a a null-measure set such that the resulting topology is trivial. This 
can be done by removing the set H (red zones in the right figure), resulting in a star-shaped open set in R 3 . 


Sublemma 1 Let n* denote the restriction of the map 
n fo B-. Then for each i we have that is a diffeomor- 
phism between B' and n(B') 

Proof. Since by definition k G BJ 4 k ^ F we have 
J n (k) 7 ^ 0 for all k G B'. Since B' is connected and n is 
analytical, we have the thesis. ■ 

Sublemma 2 We have the following inclusions: 

1. n(B[) C U, 

2. dn(B') C R. 

Proof. Let us start with the proof of item[^ By explicit 
computation we have |n(k)| 2 = 1 — A 2 (k) < 1 which 
implies that the image of n is contained in U. 

We now prove item [ 2 ] Thanks to Lemma [l] we have 
that n(B') is open. On the other hand, since n is con¬ 
tinuous and B' is compact, we have that n(B') is com¬ 
pact and then it is cl osed. Then_the trivial inclusion 
n(B') C n(B') implies n(B') C n(B'). By definition we 
have n(B-) = n(B-) U <9n(B-) with n(B') n Sn(B') = 0 
and n(B') = n(B-) Un(SB'). Then the inclusion n(B') C 
n(B') implies Sn(B') C n(3B'). Since SB' C F we have 
n(SB') C n(F). One can then verify by direct computa¬ 
tion that n(F) C R thus proving the thesis. ■ 

We now recall a result of basic topology which will be 
useful in the following. 

Sublemma 3 Let A and B be open sets such that A C 
B . Then there exists a point p such that p G inti? and 
p ^ A. 


Proof. Let us suppose that B C A. Since B is open 
and A is closed, we have B C A which contradicts the 
hypothesis. ■ 

The following result will be of crucial importance. 
Sublemma 4 The following identity holds: 

n(B') = U. (C 2 ) 


Proof. First we prove the easiest inclusion n(B') C 
U. From item [l] of Sublemma [ 2 ] we have that n(B') C 
n(B') C U (the first inclusion is trivial). Reminding that 
n(B') is open we have n(B') C U. 

We now prove that U C n(B'). By contradiction, let us 
suppose that the strict inclusion n(BJ) C U holds. Then, 
thanks to Sublemma[3j we findp G U such that p qL n(B'). 
Moreove r we c an find an open neighborhood N of p such 
that Nfln(B') = 0 and then without loss of generality we 
can suppose that p G Qk Since R has no interior points, 
n ( B i) cannot be included in R, whence n(B') D Q' is not 
empty. Let us now fix a point q G n(B') D Q'. Then, for 
any continuous path 7 connecting p and q there exist t' 
such that y(t') G <9n(B'). From item [2] of Sublemma [2] 
we have 7 (t f ) G R. Since this conclusion contradicts the 
fact that Q' is connected, we have proved the thesis. ■ 

As a consequence we have 

Corollary 1 The following inclusion holds Q' C n(B'). 

Proof. From Lemma[4]we have Q'UR = n(B') U9n(B'). 
Reminding that Q' fi R = 0 = n(B') D 9n(BG and the 
inclusion <9n(B') C R, proved in Sublemma [21 we have 
the thesis. ■ 
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We now turn our attention to the regions B'. Our first 
objective is to determine how many different B- regions 
are. The answer is provided by the following result. 

Sublemma 5 The regions B' are in one-to-one corre¬ 
spondence with the solution of the equation | A(k) | 2 = 1. 


Proof. We proved that the map rp defines a diffeomor¬ 
phism between B' and the set n(B') C P which includes 
the origin. Then, for each B' there exist a point k such 
that n(k) = 0 and it is unique. Since n(k) = 0 if and 
only if |n(k)| 2 = | A(k) | 2 — 1 = 0 we have the thesis. ■ 
Thanks to this result it is sufficient to find the solutions 
of | A(k) | 2 = 1 in the Brillouin zone. One can easily check 
that there are only 4 solutions and then 4 different regions 


We can now prove Eq. (B7) of Lemma [T] 


Sublemma 6 The region B' are given by 


Bq := {k G B|A(k) > Q,cos{2k y ) > 0}, 
:= {k G B|A(k) < 0, cos(2 k y ) > 0}, 
B' 2 := {k G B|A(k) > 0, cos(2 k y ) < 0}, 
B3 := {k G B|A(k) < 0, cos(2 k y ) < 0}. 


(C3) 


Proof. Let us denote with B^ the regions defined by 
the right hand sides of Eq. |C3| One can immediately 
see that: i) the B^ are open sets, ii) the Eb are mutually 
disjoint and iii) the union of the B^ is the union of the 
We now prove that for all B^ there exist a unique B' 
such that B^ C B'-. This fact, together with the previous 
properties of the B^ gives B^ = B • 

Clearly for all B^ we must have B^ C B'^ U • • • U B'- fc for 

some k > 1 Let us suppose then that there exist B^ such 
that B^ C B' x U • • • U B'- fe with k strictly greater than 1. 

Since we have as many B^ as B-, there must exist a B- and 
two points k a G B a and k^ G B a such that k a ,kfr G B'. 
Since B • is connected there must exist a path connecting 
k a and k& that entirely lies within B-. On the other hand, 
since the B^ are disjoint, this path would cross the border 
of the B a but this contradicts the fact that the border of 
the B a are not included in B-. ■ 

Finally, we can give the complete characterization of 
the sets n(B'). From Corollary [l] we have the inclusion 
Q' C n(B'). Since |n(k)| = 1 A(k) = 0, we 

know that the S 2 (j£_ n(B'). It is easy to check that also 

the points p± := ( 0 , ±^, 0 ) are not included in the set 
n(B'). For any region n(B'), we will determine which 
ones of the 8 open arches defined as 

E f := e ± (Lj ) 

Li '■= (0, f) L 2 :={ §,tt) (C4 ) 

L 3 :=(- f, 0 ) i 4 :=(- 


Let us consider the sets E+. If we for some t we have 
n(k) = e_|_(t) and cos( 2 k y ) 7^ 0 then it must be 


k x — f + 


k z = k x + rmr 


(C5) 


for n and m integers. Eq. (C5) then implies 


A(k) = (-l) m -(cos(fc y ) - sm(k v )). 


(C6) 


From Eq. (C 6 ) we have 


A(k) > 0 


m even A - \r: < k y < 
m odd A < k y < |7r 


(C7) 


Then, if we assume k G Bq we must have 

m even ,-^ 7 r < k y < ^tt 
m odd , §7r < k y < j7T 


(C 8 ) 


However, since the two sets of k are related by a transla¬ 
tion of (0, hr, hr), l G Z they actually represent the same 
set in the Brillouin zon e. So it suffices to consider just 
the first set in Eq. (C 8 ), that can be written as 


kx 

k z 


\ir + n\ 7 T 


\ir < k y < \ 7 r 


\ir + n\ 7 r + rmr = \ir + n\ir 


(C9) 


where we used Eq. (C5) and in the second equality in 
the second line of Eq. (C9) we used the hypothesis that 


m is even. Using again the fact that we identify points 
related by a translation of {hr, 0, hr), l G Z we find just 
two inequivalent sets 


Zi := 


k x = J7T 

\lT <ky < \tt 

k z = 


Z2 := 


kx = -jTT 

\-k <ky <\tt . 

k z = -|tt 


It is now easy to show that the images of this two sets 
under the map n are n(Zi) = Ej and n(Z 2 ) = E4 . By 
applying an analogous line of reasoning one can prove all 
the following inclusions 


(E+ U E+ U Er U E 3 ) C n(B') 
(E+ UE+U E 2 - U E 4 -) £ n(B') 
(E+U E+U Er U Eg ) c n(B') 
(E+UE+UE 2 -UE 4 -)^n(B') 
(E+ U E 3 + U E 2 - U E 4 -) C n(Bi) 
(E+ UE+U Er U Eg") £ n(Bi) 
(E+ U E 3 + U E 2 - U E 4 -) c n(Bg) 
(E+ U E+ U Er U Eg") <t n(Bg). 


(CIO) 


This result completes the proof of Eq. ( |B 8 | ) of Lemma [l] 


are included in n(B'). 
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Appendix D: The function g(cj,m) 


In this section we now show how it is possible to define 
a function m) such that the map Af defines a dif- 
feomorphism between N and the null mass-shell To- The 
purpose of the following construction is to reduce the re¬ 
gion N to a star-shaped region N by removing a null mea¬ 
sure region, and to define the function g(u, m) in such 
a way that the map Af is invertible on N. Since multi¬ 
plication by g(w, m) rescales the four vector (sincj,m) 
without affecting its direction, in order to have an in¬ 
vertible map Af it is sufficient to ensure that g(u,m) is 
radially monotonic versus m, namely g(uj,rmo) must be 
monotonic versus r. 

Let us d enot e E± the ellipses defined by the parametric 
equations ( |B5| ). We define the polynomials 

hu{r,6,4>) :=1 - r 1 2 , 

h E (r, 9,4>) :=(cos 2 * (</>) - sin 2 (0)) 2 + (Dl) 

+ (| - r 2 ( 1 - cos 2 (9) sin 2 (4>))) 2 , 

where we used the spherical coordinates m x = 
r cos#cos0, m y = rsin#, m z = r cos#sin</>. Clearly 
we have h\j(m ), > 0 Vn E, = 0 m E S 2 , 

/iE(m) = 0 m E E + U E_ and are analytic on 

U \ H. Since U \ H is star-shaped we can define 


£l(r, 9,4>) := r 


ha / 

Jo VM S , 


+ 


9,cf)) h E (s, 9,4>) 


(D2) 


The condition h\j(m ), h e(ih) > 0 Vn E U\H implies that 
the function g(r, #,</>) is radially monotonic on U \ H. 


Since h\j(m) are /ie(iti) are analytic on U \ H we have 
that g(r, #, (f) is analytic on (U \ H) \ 0. Moreover, since 
g(r, #, 0 ) is even in r we have that g(m) is analytic on 
the whole domain U \ H. Finally it is easy to check that 
g( m) diverges to +00 as m approaches the boundary of 
U \ H, and that (Vg)( 0 ) = 0. Let us define 

g(m) := g(m) + 1. (D3) 

We now check that, with this definition of the map g(m), 
the map Af defines an analytic diffeomorphism between 
U \ H and To with the property Jjg{ 0) = /. Clearly M 
is analytic in N so we just need to prove that it gives a 
bijection between N and r 0 . Let us fix a versor j in M 3 . 
Then, in the j direction we have 


•V(w,rj) = g(r,9 j,<£j) 


sin(cj) 

r j 


(D4) 


Since g(r, #j, 0j) is monotone we clearly have that A f is 
injective. We now prove the surjectivity of A/" / (cj,r) := 
W(cj,rj) on the set K := {(po,pi) E M 2 s.t. Po — Pi = 0}* 
Let us fix a point p = (po,Pi) ^ K. Since p(r, #j,0j) is 
monotone and surjective on [ 1 , + 00 ) we can find a value 
r p such that r p g(r pi #j, 0j) = p\. Clearly, since \r p \ < 1, 
also the equation sin 2 (u;p) = \r p \ 2 can be solved and then 
J\f'(uj p: r p ) = (pojPi)- Since the surjectivity of Af holds 
for any direction j, we have that Af is a diffeomorphism 
between N and To. Finally, since g(0) = 1 and Vg( 0 ) = 
0 , we have that the Jacobian of the map Af is the identity, 
which proves Eq. (B 12 ). Obviously, our choice of the 


map g( m) is not unique. 
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then h E S_. Then, let L be a group automorphisms of 
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h es h es 






